AD-A0S3  573  WISCONSIN  UNIV-MADISON  THEORETICAL  CHEMISTRY  INST  F/6  20/8 

SEMICLASSICAL  THEORY  OF  ROTATIONAL  AND  VIBRATIONAL  EXCITATION  I— ETC(U) 
JUL  77  R D OLMSTED*  C F CURTISS  AFOSR-75-2850 

UNCLASSIFIED  WIS-TCI-574  AFOSR-TR-78-0764  NL 


END 

DATE 

FILMED 


6-78 


DDC 


/ JEMICLASSICAL  THEORY  OF  ROTATIONAL  AND^IBRATIONAL 
RXCITATION  IN  COLLISIONS  OF  ATOMS 
WITH  RIATOMIC  MOLECULES^ 


X ( ] R.  D./Olmsted^-aad  C.  F. /Curtiss 


fj 


Theoretical  Chemistry  Institute  and  Department  of  Chemistry 
University  of  Wisconsin,  Madison,  Wisconsin  53706 


JJ4\  WI 
\ 


JU  > 

7 • JT  • I 


t c c '■ 


7V 


ABSTRACT 


£0  w 7 

I /;  V / 


A theory  of  rotationally  and  vibrational ly  inelastic  collisions  in 
atom-diatomic  molecule  scattering  is  presented  in  which  both  rotations 
and  vibrations  are  treated  semiclassically.  The  formalism,  however, 
allows  for  an  unambiguous  specification  of  the  initial  and  final  quantum 
states  even  though  the  S-matrix  itself  is  calculated  in  the  classical 
limit.  This  development  is  particularly  suited  to  the  description  of 
inelastic  collision  processes  among  the  highly  excited  rotational- 
vibrational  states  which  lie  near  the  continuum. 
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In  the  generalized  phase  shift,  GPS,  treatment  of  rotationally 
inelastic  atom-diatom  collisions,  the  S-matrix  is  transformed  (through 
a transformation  involving  vector  coupling  coefficients  and  representation 
coefficients)  in  such  a manner  that  the  indices  describing  the  final 
quantum  states  of  the  rotor  are  replaced  by  a set  of  Euler  angles  which 
in  the  semiclassical  regime  become  continuous.  Here  we  present  an 
analogous  treatment  of  the  coupled  rotational-vibrational  problem.  The 
S-matrix  which  is  labelled  by  discrete  indices  describing  the  initial 
and  final  vibrational  and  rotational  states  of  the  diatomic  molecule  is 
transformed  (through  a similar  transformation  which  now  involves,  in 
addition,  the  vibrational  wave  functions  of  the  isolated  diatomic 
molecule)  to  one  in  which  the  indices  associated  with  the  final  states 
are  replaced  by  three  Euler  angles,  as  in  the  GPS  development,  and  a 
continuous  parameter  associated  with  the  vibrational  motion.  This 
parameter  is  the  internuclear  separation  coordinate  of  the  diatomic 
molecule.  The  resulting  equations  are  especially  suited  to  a 
semiclassical  analysis  paralleling  that  given  earlier1  for  rotationally 
inelastic  scattering.  Since  the  present  development  parallels  that  in 
reference  1,  equations  of  that  paper  are  denoted  by  a prefix  I,  i.e. 

Eq.  (I-  ). 

This  semiclassical  development  is  believed  to  be  well  suited  to  the 

description  of  vibrational  transitions  among  the  highly  excited  states 

of  the  diatomic  molecule.  In  particular  it  is  applicable  to  such 

2 

problems  as  the  termolecular  recombination  of  hydrogen  via  the 
mechanism 
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H + H --  H* 

H*  + M - H2  + M . 

The  stabilization  of  the  metastable  by  collision  with  the  third  body 
M (M  = Ar,  He,  etc.)  involves  both  vibrational  and  rotational 
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de-excitations  within  the  levels  close  to  the  continuum.  The 
application  of  this  semiclassical  method  to  this  problem  is  currently 
being  investigated. 

The  problem  considered  is  the  nonreactive  collision  of  an  atom  and 
a diatomic  molecule  interacting  through  a Born-Oppenheimer  potential. 

A center  of  mass  coordinate  system  is  chosen  in  which  § is  the  vector 
between  the  atoms  comprising  the  molecule  and  r is  the  vector  between 
the  atom  and  the  center  of  mass  of  the  diatomic  molecule.  The  internal 
state  of  the  diatomic  molecule  is  described  by  the  wave  function 

Z(v*;0  y’J(E) 


and  is  labelled  by  the  vibrational  quantum  number  v and  rotational 
quantum  numbers  l and  m . The  Y™(£)  are  normalized  spherical 
harmonics  and  the  radial  functions  Z(v£;£)  are  eigenfunctions  of 
the  differential  equation 


r *£  d2 
L"  2M  d^ 


flilU+Ii 

2MC2 


+ V.  + 


- Eb(v£)]Z(vJt;0  = 0 . 


(1) 


w 
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In  this  equation  M and  , respectively,  are  the  reduced  mass  and 
potential  energy  of  the  isolated  diatomic  molecule  and  Eb  (Eb  < 0 ; 
measured  from  the  separated  atom  limit)  is  the  binding  energy  of  the 
molecule.  These  radial  functions  satisfy  the  completeness  and 
orthogonality  conditions  that 

l Z(v£;£)*  Z(vt;e')  = 6(5-?') 
v 

and 

| Z(v*;C)*  Z(V£s5)d5  = . 

An  integration  over  continuum  states  is  implicitly  contained  in  the 
completeness  relation  Eq.  (2). 

All  of  the  information  about  the  scattering  process  is  contained 
in  a set  of  complex  functions,  Q^(vfcX;S£r)  , which  satisfy  the  following 

4 

pair  of  coupled  equations; 


(2) 

(3) 


~ exp[iQ(  + )(vU;SEr)]  = ^ 


dS’  d?'  dS"  dC 


F(±)(n;SES'C'  ;r)  V(1  ^rC'0') 

(4) 

{F(+)(n;S"E"SlCl  ;r)*  exp[iQ(  + )(v$X;S"Cr)] 

- F(‘)(n;S'T,S,r;r)*  exp[i(/"^  (vlX;S"5"r)]}  . 


This  pair  of  equations  is  obtained  by  suitably  transforming  those  obtained 
from  an  integral  equation  based  upon  a separation  of  the  total  atom- 
molecule  interaction  potential  V(r£0)  into  elastic,  V^(r)  , and 
inelastic,  V^(rf;0)  , contributions.  The  angle  0 is  that  between  the 


4 

vectors  r and  § . The  matrix  elements  F^^fcAiSCS'C  ;r)  are 
defined  by 

F(±,(jU;S£S'e';r)  = £ X(^;v*XL;SC)»  f(±,(v£A;r)  xU^AUS'C' ) 

vilAL 

(5) 

in  terms  of  the  basis  functions 

X(U;vJUL;SO  = (-1  )L  (8tt2)  '*  [(2U1 ) ( 2X+1  )(2L+1  )]/z 

(l  l L ) fX  X L]  , (6) 

l (-i)  y DLCS)  Z(v£;5)*  . 

vy  [0  v -vj  [0  -y  yj  M 

The  f^(v£A;r)  are  solutions  of  the  radial  wave  equation  involving  only 
the  elastic  part  of  the  interaction  potential.  They  are  defined  in 
detail  in  reference  4. 

The  coupled  equations,  Eq.  (4),  may  be  written  in  a more  compact 
form  by  defining  an  operator  F^  by  its  action  on  an  arbitrary 
function  Y (ft)  as 

F(±)  Y(n)  5 | F(±)(£A;m';r)  Y(fl' ) dfl'  , (7) 


where  fl  stands  for  the  combined  indices  S£  and  dn  = d£  dS  . It 
follows  that 


5 


Another  operator,  G , is  defined  by 


G(±)  Y(fl)  = 


;r)  Y(n')  dn'  , 


(9) 


where 


G(±)(U;flft';r)  = l x(SX;v£AL;n)*  [^r  f(±,(v£X;r)]  x^vULifl' ) 


vJtXL 


(10) 


It  follows  immediately  from  these  definitions  that 


,(±)  = F(±)-l  F(±)'  = FU)'  F(±)-l 


and 


F(±)  G<*>  - G<+-) 


(11) 


(12) 


where  a "prime"  indicates  differentiation  with  respect  to  r . The 
equations  for  and  can  be  uncoupled  to  give  the  second 

5 

order  differential  equations 


W Q 


- (ft  £ Q)2  = 2y  e"iQ  F V(1)  F _1  eiQ 


3r 


+ ft2  e“iQ[G+  F(^r  In  V(1})  F _1 


(13) 


+ F V(1)  G V(1,_1  F’1]^  elQ)  , 


where  (to  simplify  the  notation)  the  superscipts  (±)  have  been  omitted. 
This  expression  may  also  be  written  as 


ifi2  Q - (ft  £ Q)2  - 2y  e'iQ  l / eiQ 


+ ift2  e'lQ[G  + U if1  + V G l/_1]elQ  Q , 


(14) 
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where  U and  1/  are  two  additional  operators  defined  by 


and 


V = F V^1  V"1 


U = F(^rV(1))F_1 


(15) 

(16) 


Useful  properties  of  the  operators  F and  G are  obtained  in  a 
manner  analogous  to  that  discussed  in  reference  1.  To  this  end, 
operators  K and  J are  defined  as 

K 


and 

J 

The  angular  momentum  operators,  L.  and  M.  , are  those  used  in  I.  The 

J J 

factors  K and  J are  given  in  terms  of  the  initial  angular  momentum 
quantum  numbers  as 


UM52)'1  U,2*i22-t32-2l  K(^)~n  L3"  Lz  L,-] 

1,2  d2  " <>7) 


*VE 


2M  dCz 


= M,2  + M22-«32+2l  |^lij(W-"«3n«2«3n.  (18) 


K = 1(1+  l)-h2 

and 

J = X( A+l  )fi2  , 


(19) 

(20) 


while  l is  defined  in  terms  of  the  total  initial  binding  energy  of 
the  diatomic  molecule,  (Eb  < 0 ) , by 


■ 
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E = Eb(^}  - w1  *a+1)  • (21) 

These  K and  J operators  have  the  very  interesting  and  useful 
properties  that  x is  an  eigenfunction  of  K * and  J*  with  the 
eigenvalues  being  the  energy  inelasticity  and  the  change  in  the  relative 
angular  momentum,  viz. 

K*  xliX;v*AL;SO  = -(E-E)  x(AA;v£AL;SO  (22) 

and 

J*  x(*A;v£AL;S£)  = h2[A(A+l ) - A(A+1 )]  x(£A;v£AL;S£)  . (23) 

The  energies 

fi2fc2 

E ’ TT  ' <24> 

and 

E = Et  - Eb(v£)  (25) 

are  the  relative  translational  energies  of  the  molecules  before  and 

after  a collision;  E^  is  the  total  center  of  mass  energy. 

( + ) 

The  functions  f'-'(v£A)  central  to  the  operators  F and  G can 
be  written  in  the  form  (see  Eq.  (1-39)) 

f(±,(vU)  = (^-)  exp[J  S(±)(v£X)]  . (26) 

The  functions  S^(vlA)  , or  dropping  the  superscripts,  S(v£A)  , 
depend  on  v,  £ and  A only  through  the  quantities  E and  A(A+1  )"h2  . 

iL _ _ 
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Therefore,  they  may  be  expanded  in  a double  Taylor  series  about  their 
values  at  E = E and  A(A+l)fi2  = J ; 


S(vJX)  ' l ij7rSmn(E-E>m  [MA*m2-J]n 


where 


~m  ~n 

S = 2—  S(v£A)  . 

mn  9Em  9Jn 


These  properties  enable  the  integral  operators  F and  G to  be  written 
in  differential  forms.  By  the  definition  Eq.  (8),  F acting  on  an 
arbitrary  function  Y(&)  is  given  by 

F Y(fl)  = l xU*»vUL;R)*  f(v£A)  x(^;vUL;R' ) Y(ft' ) dfl1 
v£AL 

(29) 

= (&  4 l x(Z5i;vJlAL;n)*  e(l/fl)S^vU)  x(^;vUL;R' ) Y(R')  dft1  . 
v£AL 

From  the  eigenvalue  relations  Eqs.  (22)  and  (23)  and  the  expansion 
Eq.  (27),  it  follows  that 


,f>2N1/4 


F Y(fl)  = Ip”  l xUX-.vJLALifi)*  j Y(fl')  exp[i  l S|nn(K*),,,(J*)M] 


vm  / . n-i 


x(£A;v£AL;fi' )dfl'  . 


Since  the  operators  K and  J are  hermitian  and  the  x(^Jv£AL;R)  obey 
the  completeness  relation  that 


l x(ZX;v£AL;fl)*  XtfS;v£AL;fl')  * 6(n-fl')  = 6(S-S' ) 6(5-?' ) , 

vUL  (31) 
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Eq.  (30)  is  equivalent  to  the  expression 


where 


F Y(Q)  = (|y‘)  exp[^S]  Y(fi) 


s = y (--J-L-  s k™  Jn 
L m! n!  mn  h ' 


Because  Y(ft)  is  an  arbitrary  function,  F has  the  differential  form 


*•<&>*  s3 


Similarly,  6 can  be  expressed  as 


G = IS'  . 


A semiclassical  solution  to  the  equation  for  Q can  be  obtained 
after  suitable  approximations  are  made  concerning  the  operator  F . 
First,  the  function  , defined  by  Eqs.  (1-37)  and  (1-38),  is 

written  in  the  series  form 


S(0)(vU)  = si0)  + i"h  sj0)  - f)2  si0)  + ...  , 


where  the  are  given  by  Eqs.  (1-53)  through  (1-56).  In  a similar 

manner,  the  functions  S^(vU)  [or  S(vU)  ] and  the  operators  K 
and  J have  the  series  representations 


The  operator  S defined  by  Eq.  (33)  may  now  be  expanded  as 


S = l (ill)  Sn  , (47) 

n=0 


11 


s 


n 


n 


I s 

k=0 


nk  * 


(48) 


Each  operator  S ^ is  at  most  of  order  k in  derivatives  with  respect 
to  angles  and  the  internal  separation  distance  of  the  diatomic 
molecule.  Specifically,  the  first  few  terms  are 


where 


00  / i 

: = Y L~.V  . k m c 

nO  m^0  m!  *0  bn,mp 


(49) 


and 


nl 


n2 


00  i i 

L m!  K0  ^Sn-l,m,l  J1  ' Sn-l,m+l,0  Kl^  (50) 


m=0 

00  / , \m 

m=0 


00  / i \m  9 

L m!  ‘o'U  Sn-2,m+2,0  K1  + ^ Sn-2,m,2  Jl' 


" Sn-2,m+l ,1  J1  K1  + Sn-2,m,l  J2] 


(51) 


00  / t \m  m-1  , . . 

+ V \ / v i/  & ls  i/  ni-4  ~ 1 q 

mil  »<  Jo  0 2 0 


n-2,m,0 


S = (-L)m  (J_)n  s 

P.rryi  V bp 


(52) 


Except  for  the  last  term  in  Sn2  , the  infinite  sums  that  appear  in  the 
Snk  have  the  effect  of  just  replacing  the  initial  translational  energy 
E , which  parameterizes  the  Sp  mn  , by  the  effective  energy 

Et"W"V(^  ; 
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® / , xtn 


(-1)  K m 

m=0  m!  0 Jn,n**,t  m=Q 


y a.-ji  K s = y 

£ _ ml  ^n_m+A  t l*  - HI ! 


/ i \in 


j0  iT»t  - & - ««  - ^ <^>"  Sn,i(t  <«) 


= Sn  At  I - K = SU)_  • 

^ ! E=Et  ■ 2Hf  ' v(c)  ^ 

In  tenns  of  the  S(£)n  just  defined,  the  Sn(<  are  given  by 


Sn0  = S(?)n,00  S^n  ’ 


(54) 


Snl  = "S(C)n-l,l,0  K1  + S(C)n-l,0,l  J1 


(55) 


Sn2  = S^n-2,C0  J2  + ^S^n-2,2,0  K1  " S(?)n-2,1)1  K1  J1 

+ 1/2S(S)n-2,q2  J]Z  ' S(^n-2,1}0  KZ  " ^ 

" ^6S(On_2>3p[Ko’^o,K2^  ’ 


where  the  results  of  the  Appendix  have  been  used  in  the  expression  for 
Sn2  . The  expressions  for  SnQ  and  Snl  are  identical  in  form  to 
Eqs.  (1-68)  and  (1-69)  with  Sn  4t  replaced  by  the  corresponding 
C-parametrized  quantities  SU)n  . The  expression  for  Sn2  contains 
additional  terms  involving  commutators  of  Kg  with  K2  . 

In  this  way,  the  central  quantities  that  must  be  calculated  for  the 
description  of  the  full  vibration-rotation  inelastic  scattering  problem 
are  reduced  to  expressions  similar  to  those  arising  in  the  rotation 
problem  alone.  The  difference  is  that  now  the  generalized  action 


13 


integrals  are  parameterized  by  the  internuclear  separation  distance  of 
the  diatomic  molecule.  This  only  affects  the  effective  translational 
energy  at  which  the  generalized  action  integrals  are  evaluated. 

Using  the  expansion,  Eq.  (47),  one  finds  that 

e(i/fl,S  = H (ifi)n  Rn]  e'Sn  e"Sl  e(i/f,)S° 

n 

and 

e-(i/fi)S  = e-(i/fi)S0  eSj  eSnQ  (ifl)n  _ 

n 

The  operators  Rn  and  Pn  are  given  by  Eqs.  (1-75)  through  (1-81),  but 
once  again  the  appropriate  ^-parameterized  quantities  must  be  used. 

The  problem  is  thus  reduced  to  that  of  reference  I and  the  further 

analysis  is  not  repeated  here.  The  essential  difference  is  that 

the  various  quantities  have  an  additional  dependence  on  the  variable  £ . 
The  result  is  that  when  Q is  written  in  the  semiclassical  form 


(57) 

(58) 


q ■ ij  mnQn . 


(59) 


the  coefficients  Qq  and  Q1  satisfy  the  equations 


<F  V2  - 600<£  V + 2“  *00  ■ 0 


(60) 


and 


2t®00-<£  Wfqi>  = -2u  V^Qo 


+ {Goo  woo  + woo^Girl/oo^  + ti,oo*-l/irGoo-*"  ZGi^ar  V 


(61) 
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The  expressions  for  the  quantities  6,  P and  W are  the  same  as  those  in  I. 
It  is  remarked  once  again  that  all  of  these  quantities  depend  on  C 
through  the  S(£)  .*  which  are  central  to  their  definitions.  It 

11)^  l 

should  also  be  noted  that  implicit  in  the  expressions  for  G, 

1/  and  W is  still  given  by 


V(1)  = V(1)(S',C)5  (62) 


but  the  angles  and  0 used  to  specify  the  rotation  S ' , 


L 


S'  = S(tt,0,tt)  , 


(63) 


depend  on  C as  well; 


<P  = = - 


K 


V2 


MC  S(C)0,1,0  ’ 


0 = 0(C)  = -2jy2  S(C)QjQ>1  . 


(64) 

(65) 


To  first  order  in  the  inelasticity,  the  solutions0  of  Eq.  (60)  are 
given  by 


Q^lvUSCr) 


In  Z(v£;C)  + 2n(vl)  f° |f^(vlA)|2  V*1  ^T^.Odr' 

' r 

(66) 


and 


q£+)(vSiASCr) 


= In  Z(v£;C)  + 2n(vfc)  - & f |f(_)(v£A)  |2  V(1  ^T^.Cjdr' 

* JrnU) 


£r  f"  |f(+)(vH)|2  vn)(T(+),C)dr' 


r0(S) 


(67) 


- — - 
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The  rotations 

T^  are  given  by 

T^  = (tt,^(^),tt)  S(tt,0(C),it)  , 

(68) 

and 

T<->  - 

(tt,2tP(0)(C)  - iMS)  ,tt)  S(7t,29(0,U)  - 8(5)  ,ir)  , 

(69) 

where  i(/0)(C)  and 

9^(0  are  the  values  of  MC)  and  9(C) 

at  the 

turning  point  r^C)  . The  results  Eqs.  (66)  and  (67)  may  then  be  used 
in  Eq.  (61)  to  generate  the  solutions  and  . 

Knowledge  of  the  Q^(vUSO  is  sufficient  to  determine  the 
S-matrix  (in  an  unusual  representation) 

S(v£XSC)  = lim  eiq(  f (70) 

f -►  » 


in  which  the  initial  state  of  the  system  is  specified  but  the  final 
state  indices  are  replaced  by  the  angles  S and  the  distance  C . 
From  this  S-matrix,  the  transition  amplitude,  T(v£A;vJlAL)  , may 
readily  be  determined; 


5ij  o “ V*2 

T(vlX;vUL)  = 6(v£A0| v£AL)  - (-1)  A (8tt  ) 


j x(^X;v£XL;SC)  S(v£ASC)dC  dS  . 


(71) 


These  expressions,  Eqs.  (70)  and  (71),  for  the  S-  and  T-matrices  are 
exact  but  when  combined  with  the  previous  results  for  Q(v£,XSE) 
provide  an  approximation  in  which  the  rotational-vibrational  scattering 
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is  treated  in  a semiclassical  manner.  Despite  the  semiclassical  nature 
of  the  development,  the  initial  and  final  quantum  states  are 
unambiguously  specified  in  the  scattering  amplitude. 
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APPENDIX 


Since  Kg  and  K ^ do  not  comnute,  the  sums  in  the  expression 


y (-1  )m  my1  K a K ^m-A-1  . 

ii  ""  4=o  0 20  "-2-m0 


(which  is  the  last  term  In  Sn2,  Eq.  (51))  cannot  be  evaluated  as  in 
Eq.  (53).  However,  these  sums  may  be  evaluated  in  the  following  way;y 


V V (~1 )"'  k 6k  j^ni-a-1  <.  y y (-1  )t+^  K 6 K *r^-4  c 

Jo  m!  *0*2*0  bn-2,m0  " t=0  Tt+TTT  0 2 *0  bn-2,t+l,0 


? f i-i 

4=0  p=o 


p'.a!  -fr  K9  ttt 


4!  2 p!  n-2,p+4+ip 


Y d!4!  ('K0  u,  ("K0  aT,P  c 

'/p  Tp+4+tyi  4!  K2  p!  n-2,10 


fl  . p ("Kg  -j^)4  ('Kg  ^)P 

-I  da(l.a)AaP-^<2-^L_Sn.2J0 
o 


.’l  e^"1 


-1  ;K0(3/3E)  . -aK0O/9E) 
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where  the  substitution  p = t-4  , the  integral  representation  for  the 

O 

beta  function  p!4!/(p+4+l ) ! , and  the  definition  of  the  S„  , 

r r p,mn 

Eq.  (52),  have  been  used.  The  integral  is  evaluated  by  defining  a 


comutator  superoperator  A as 
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